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Abstract 

Interface technology for geometrically nonlinear 
analysis is presented and demonstrated. This tech- 
nology is based on an interface element which 
makes use of a hybrid variational formulation to pro- 
vide for compatibility between independently mod- 
eled connected subdomains. The interface element 
developed herein extends previous work to include 
geometric nonlinearity and to use standard linear and 
nonlinear solution procedures. Several benchmark 
nonlinear applications of the interface technology are 
presented and aspects of the implementation are dis- 
cussed. 

Introduction 

Rapid detail modeling and reliable analysis 
methodology are needed in order to provide tractable 
computational capabilities early in the design of 
aerospace and ground vehicle structures. Although 
critical structural details are often not well-defined in 
the early preliminary design phases, achieving an 
understanding of the effect of certain details can 
greatly impact early design decisions. In addition, 
rapid modeling and analysis capabilities that include 
such design detail provide early insight into and pre- 
dictions of the failure mechanisms of structural con- 
cepts. Design verification, often satisfied through 
refined analysis, may be impacted significantly by 
computational methods, simulations and techniques 
which predict accurately detailed stress states such 
as yield and fracture stresses. Detail stress analysis 
often relies heavily on the finite element method 
which requires fine mesh discretizations of the struc- 
ture in the vicinity of stress or strain concentrations. 
The resulting finite element models generally con- 
tain both fine and coarse mesh discretizations. Tedi- 
ous, costly and time consuming transition modeling 
is generally used between the subdomains of fine and 
coarse refinement. Such transition modeling is also 
typically used for assembling independent substruc- 
ture models, inserting part or component detail, 
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global/local modeling and modifying structural fea- 
tures. The extensive labor required for such transi- 
tion modeling often precludes the use of the finite 
element method early in design. 

Recently, a method for connecting finite ele- 
ment models without the use of transition modeling 
has been developed. 1,2 This method, called interface 
technology, is an improved technique for connecting 
multiple dissimilar meshed subdomains or substruc- 
tures to form a single finite element model. Com- 
pared with conventional transition modeling tech- 
niques, this method allows designers to create mod- 
els for structural analysis faster and with less effort 
while retaining computational accuracy. The method 
promises to significantly reduce engineer modeling 
time and effort, thus enabling comprehensive early 
design. 

Unlike many global/local 3,4 and substructuring 5 7 
approaches, the interface technology provides a 
means of connecting independently modeled sub- 
structures or subdomains whose nodes along a com- 
mon boundary need not coincide. The technology is 
based on a variational formulation which enforces 
the compatibility between the connected subdomains 
in a weak sense. This compatibility provides for the 
interaction between the independently modeled sub- 
domains. Conventional global/local techniques, in 
which global coarse models are used to provide 
loads or displacements 8 for local refined models, do 
not account for full subdomain interaction including 
changes in the load path prediction due to the refined 
local modeling. In addition to global/local modeling, 
interface technology enables assembly of independ- 
ently created model components, repair patch model- 
ing, contact/friction/sliding applications and detail 
model insertion. 

The interface technology has been shown to 
maintain solution accuracy for a wide range of appli- 
cations. 1,2,9,10 The method was developed and vali- 
dated on benchmark problems in reference 1, refor- 
mulated in the form of an element in reference 2, 
and demonstrated on large scale structural applica- 
tions in references 9 and 10. The interface technol- 
ogy has been incorporated into the commercially 
available finite element computer code, 
MSC/NASTRAN. 11,12 However, until recently, the 
interface technology has been limited to linear stress 
analysis. Although linear analysis is quite useful, it 
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can miss important phenomena. Ignoring nonlinear 
effects can lead to nonconservative predictions for 
softening systems and severe weight penalties for 
stiffening systems. 

Thus, the purpose of this paper is to present the 
extended capabilities of the interface technology to 
geometrically nonlinear analysis. This nonlinear 
capability is developed, described and demonstrated 
on representative structures. A new interface ele- 
ment solution strategy for both linear and nonlinear 
analysis is developed and described. The effective- 
ness of the approach is demonstrated using five elas- 
tic plate and shell applications including: a pressure- 
loaded clamped square plate, a pressure-loaded 
clamped cylindrical panel, a point-loaded hinged 
cylindrical panel, a moment-loaded cantilevered 
plate and a compression-loaded composite panel. 
These application problems exhibit nonlinear re- 
sponses that verify the capabilities of the technology; 
namely, (1) softening, (2) stiffening, (3) snap- 
through, (4) snap-back, (5) large rotation behavior 
and (6) postbuckling. The guidelines for developing 
this technology include the requirements that it be 
compatible with general-purpose finite element 
codes, valid for a wide range of finite elements, cost- 
effective and accurate for predicting the linear and 
nonlinear responses as well as the overall and de- 
tailed stress states of structural components. 

Ove rv iew o f t he Interface Element 

The interface element was developed in detail in 
reference 2 and is briefly described herein. It allows 
the independent modeling of different substructures 
or components without concern for one-to-one nodal 
coincidence between the finite element models. 
Moreover, it acts as “glue” between independent 
finite element models with different mesh densities 
and nodal layouts. Thus, the interface element pro- 
vides a transition modeling role through an analytical 
variational procedure and avoids the use of transition 
meshes. It alone does not improve the quality of the 
finite element results obtained by a particular model 
but rather improves the efficiency of the modeling 
and makes effective use of existing finite elements. 
By eliminating or reducing transition modeling, the 
introduction of distorted elements, which may de- 
grade the quality of the solution, is limited to what is 
required to represent the geometry of the structure. 
In addition, the elimination of unnecessary element 
distortion errors allows the use of a less refined mesh 
and therefore, the same qualitative results may be 
obtained with fewer degrees of freedom. 

Consider, for example, the domain shown in 
Figure 1 and modeled as three independently dis- 
cretized substructures. The interface element is dis- 


cretized with a mesh of evenly-spaced pseudo-nodes 
(open circles in the figure) which need not be coin- 
cident with any of the interface nodes (filled circles 
in the figure) of any of the substructures. The hybrid 
variational formulation 1,2 employs an integral form 
for the compatibility between the interface line ele- 
ment and the finite element substructures. The dis- 
placement vector, v, of the interface element is as- 
sumed to be independent of the displacement vec- 
tors, k, of the substructures to which it is attached. 



Figure 1. Typical Interface Element 
Definition. 


Nonlinear Formulation 

In the previous work presented in reference 2, 
the interface element was limited to linear analyses. 
In the present work, the range of applicability is ex- 
tended to the geometrically nonlinear regime. The 
approximation in the domain, Q*, is made in 
terms of the displacements, u V} and v, , and the inter- 
face tractions, X iy Tractions and displacements are 

prescribed on boundaries and respectively 
(see Figure 2). 



Figure 2. Linking of Two Domains Using 
Interface Technology. 
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Using the weak form of the principle of virtual 
work, leads to 


large displacement matrix, and k T is the tangent 
stiffness matrix. 


SWfc - J ISufcT Ffc dQfc 

- j Su k T <p k dS^ - \Su k T X k dTj 

s (°) r k 

k 

+ \&k T h dr/ = o (1) 

r / 

where for the k^ substructure, F k are the external 
applied forces and <p k are the applied tractions. As 
usual, the satisfaction of prescribed displacement on 

is implied by the approximation for u { . The N ss 

subdomain equations are completed by enforcing 
compatibility through the use of the constraint inte- 
gral for the Nj interface elements 

<5^, = I / s4h-«y) dr V = o ( 2 ) 

7 =1 

where w„(i) is the number of substructures connected 
to interface element /. Use of this constraint integral 
corresponds to the “frame” method of reference 13. 
Therefore, for the entire domain 


Upon assembling the element contributions in 
each subdomain, Equation (4) may be written as 

j<j T 8edQ k =8ql(K T ) k q k 

a* 

where ( K T ) k is the tangent stiffness matrix for the 
kf h substructure. 

The independent approximations for the finite 
element displacements, interface displacements, and 
interface tractions are, respectively 

“«> = N iWj 

n = Ttfsi (7) 

% = R ij°ij 

where q tj and q Si are the nodal degrees of freedom 

corresponding to u t j and v,, and are the unknown 
coefficients of the Lagrange multipliers, ky. The 
matrix Ny is the matrix of finite element shape func- 
tions on substructure j along interface i, T t is formed 
as a result of passing a cubic spline through the 
evenly-spaced pseudo-nodes, Ry is formed as a result 
of using constant functions for linear finite elements 
and linear functions for quadratic finite elements. 


_ H ss Nj (3) 

SW = ISW* + I<5W C| - =0 
k=l i=l 

Discretization of displacements in each domain 
and displacements and Lagrange multipliers on the 
interface yields the final system of equations. At the 
finite element level, the first integral term in Equa- 
tion (1) may be written as 

Jtx r &dft* = \5 e t <j dO* 

nel{k) (4) 

= I }& / crdQ e 
n-\ Qe 

where nel(k) is the number of elements in the k** 1 sub- 
structure. Equation (4) is evaluated using the non- 
linear strain displacement relations and the constitu- 
tive relation for elastic, small strain behavior and 
may be rewritten as 

\ Se T o = SqJ (*o + k a +k L )q e 

Qe (5) 

= SqJ(k T )q e 

where q, is a vector of nodal generalized displace- 
ments for the element, k 0 is the linear stiffness ma- 
trix, k a is the geometric stiffness matrix, k L is the 


Hence, Equation (3) may be rewritten as 


Nr, 


SW= 1 Stl{K T ) k q k -Stlf k 

k = 1 

Nj t f 
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and defining the interface matrices as 

Nj jRij dry and Gy = \T? Ry d Ty (,0) 


r ij 

Thus, for arbitrary q k on £2 k 

fy n ss(*) 

Z ( K T) k <lk -fk+Y *L M ij a ij =° 

1=1 J =1 


k=\ 

for arbitrary q f on r, 


Nl «**(*) 
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for arbitrary Cty on T; 
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(13) 


Thus, the resulting system of equations is given 
in matrix form as 


(14) 


' K t 0 m{ 

* 


/ 

0 0 Gj 

9s 


0 

mJ Gj 0 

a 


0 


where K Tt q and / are the assembled tangent stiffness 
matrix, displacement vector and force vector for the 
entire structure, and A/,, G,, q s and a are the assem- 
bled Afy, G tf , q^y and Cty for all interface elements. 

The assembled tangent stiffness matrix K T is a block 
diagonal matrix containing the tangent stiffness ma- 
trices ( K T ) k of each of the substructures along its 
block diagonal. The interface element “stiffness” 
matrix and vector of unknowns are given in the same 
form as that given for the linear interface element of 
reference 2. However, in the moderate- to- large strain 
regime, the integral limits for the coupling terms of 
Equation (10) are related to the reference state (the 
deformed configuration in this implementation) and 
account for straining of the interface. This straining 
has been neglected herein but does not adversely 
affect the interface element performance for this 
study. The system of equations in Equation (14) is 
symmetric but not banded nor positive definite. 
Therefore, standard Cholesky solvers may not be 
used, unless row or column pivoting is performed to 
obtain the solution. The number of additional de- 
grees of freedom associated with the interface ele- 
ment is generally small in comparison with the total 
number of degrees of freedom in the structure. Thus, 
modeling flexibility is provided at a relatively small 
computational expense. The computational expense 
may be reduced additionally as the efficiency of new 
solution algorithms for the system of equations in 
Equation (14) is increased. This nonlinear interface 


element was implemented within a general-purpose, 
finite element code COMET/AR. 14 This implementa- 
tion is briefly described in the next subsection. 

Nonlinear Element Implementation 

The nonlinear solution procedure employed 
herein is based on a Newton/Raphson incremental 
strategy for automatic load step control. The so- 
called modified Newton/Raphson method, which 
forms and factors the tangent stiffness matrix peri- 
odically rather than at every nonlinear iteration, has 
been used in the nonlinear analyses of the bench- 
mark applications presented in the next section. A 
corotational formulation 15,16 which identifies the ref- 
erence state as the current de formed configuration is 
used to describe the motion. This formulation sepa- 
rates the rigid body motion from the strain-producing 
motion thus allowing for either linear or nonlinear 
strain-displacement relations at the finite element 
level. Deformations are computed based on the 
original configuration within the local corotated 
frame. An arc-length control strategy, 17 the most 
general approach for obtaining the load-deflection 
response, is used to handle limit points with both 
snap-though and snap-back behaviors. 

The standard nonlinear solution strategy has 
been adapted to incorporate the interface element. 
Interface element “stiffness” matrices are computed 
and assembled along with the finite element tangent 
stiffness matrices. The displacement solutions are 
used to compute the increment in the arc-length and 
predict the solution at the next load step. In analyses 
with the interface element, these displacement solu- 
tions contain coefficients of the Lagrange multipliers 
as well as the nodal displacements (see Equation 
(14)). These Lagrange multiplier coefficients are 
eliminated when computing the arc-length increment. 
In addition, in the standard procedure, the number of 
negative roots is used to determine the structural 
instability during loading as well as to identify 
changes in the loading direction. This is not a viable 
approach in this work since the indefinite system of 
Equation (14) could have both positive and negative 
roots even along the stable equilibrium path. Thus, 
an approach based on the change in incremental 
work is used herein for the interface technology 
rather than the number of negative roots to determine 
instability and change in load direction. The change 
in incremental work is essentially the so-called cur- 
rent stiffness parameter. 18 

Nu meri c al Results 

The capabilities of the interface element for 
geometrically nonlinear analysis are demonstrated on 
several benchmark applications. The application 
structures exhibit a wide variety of response charac- 
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teristics, namely, stiffening, softening, snap-through, 
snap-back, large rotation, and postbuckling re- 
sponses. These response characteristics test exten- 
sively the robustness of the nonlinear interface ele- 
ment implementation. A nine-node Assumed Natu- 
ral-coordinate Strain (ANS) shell element 19 has 
been used in each of the applications discussed 
herein. The ANS element has five degrees of free- 
dom at each node (three displacements, two rota- 
tions, and no drilling degree of freedom). Results 
from each interface element analysis (denoted cou- 
pled analysis herein) are compared with a reference 
solution obtained by using a reference model (with 
no interface) in which the discretization of the most 
refined independent substructure model extends over 
the entire domain. The reference solutions are ob- 
tained from analyses performed as part of this study. 
These solutions, in excellent agreement with those 
published in the open literature, are used to assess 
the performance of the interface technology for the 
given finite element type, mesh discretization, and 
nonlinear solution strategy. 

Clamped Square Plate 

An isotropic square plate which is clamped on 
all four edges (shown in Figure 3) and is subjected to 
a uniform pressure load is used to demonstrate the 
capability of the interface technology for a structure 
exhibiting a stiffening response. The plate material 
properties are 1000 ksi for the Young’s modulus and 
0.3 for the Poisson’s ratio. The plate edge length, L, 
is 10 in., and the thickness, t, is 0.1 in. The coupled 
finite element model is shown in Figure 4. The inter- 
face element is located at the panel midlength 
(shown by the gray shaded line in Figure 4). In 
Figure 5, the normalized pressure load is given as a 
function of the transverse deflection at the center of 
the plate (point C in Figure 3) normalized by the 
plate thickness. The response curve indicates the 
stiffening behavior of the plate as the load is in- 
creased. The coupled analysis, utilizing an interface 
element as the coupling agent, is in excellent 
agreement with the reference solution. 



Figure 3. Clamped Square Plate Subjected to 
Uniform Pressure. 



Figure 4. Finite Element Model for Clamped 
Square Plate. 



Normalized Transverse Deflection, Wf/t 

Figure 5. Central Deflection of Clamped Square 
Plate. 

Clamped Cylindrical Panel 

An isotropic cylindrical panel which is clamped 
on all four edges (shown in Figure 6) and is sub- 
jected to a uniform pressure load is used to demon- 
strate the capability of the interface technology on a 
shell structure as well as for a structure exhibiting a 
softening response. This panel configuration is given 
in the literature and attributed to Brebbia and Con- 
ner. 20 The panel material properties are 450 ksi for 
the Young’s modulus and 0.3 for the Poisson’s ratio. 
The panel planform is square with a length, L, of 20 
in. and a radius, R, of 100 in. The panel thickness, 
t, is 0.125 in. 

The finite element model shown in Figure 7 in- 
volves two independently modeled, coupled substruc- 
tures. The interface element is located along an ax- 
ial line at half the panel planform dimension. In 
Figure 8, the pressure load, p, is shown as a function 
of the transverse deflection at the center of the panel 
(point C in Figure 6). The response curve indicates 
the softening behavior of the panel as the load is 


5 

American Institute of Aeronautics and Astronautics 



increased. The coupled analysis utilizing an inter- 
face element is in excellent agreement with the ref- 
erence solution. 



Figure 6. Clamped Cylinder Subjected to 
Uniform Pressure. 



Cylindrical Panel. 



Figure 8. Central Deflection of Clamped 
Cylindrical Panel. 

As discussed in a previous section, the corota- 
tional formulation allows for the use of either linear 
or nonlinear strain-displacement relations at the fi- 
nite element level. (The corotational formulation 
includes a substantial part of the nonlinear response 


character by separating out the rigid body motion.) 
The interface technology provides the flexibility to 
use different strain approximations in each of the 
independent substructures. The effect of the nonlin- 
ear terms in the strain-displacement relations on the 
transverse deflection is shown in Figure 9. Results 
are obtained with linear and nonlinear strain- 
displacement relations and are denoted in Figure 9 
by (L) and (NL), respectively. Three reference solu- 
tions are obtained and are given by the solid lines in 
the figure. The first reference solution, denoted ref- 
erence (NL) in the figure, is obtained using nonlinear 
strain-displacement relations and the discretization 
of the most-refined substructure of the coupled model 
(see Figure 7) over the entire domain of the panel. 
This finite element model will be referred to hereaf- 
ter as the refined reference model. The second refer- 
ence solution, denoted reference (L) in the figure, is 
obtained using linear strain-displacement relations 
and the refined reference model. The third reference 
solution, denoted coarse reference (L) in the figure, 
is obtained using linear strain-displacement relations 
and the discretization of the less-refined substructure 
of the coupled model (see Figure 7) over the entire 
domain of the panel. This third reference solution is 
used to bound the results of the coupled analyses, 
and the finite element model used will be referred to 
hereafter as the coarse reference model. When lin- 
ear strain approximations are used, the predicted 
structural behavior is not as soft as that behavior ex- 
hibited with nonlinear strain approximations. This 
behavior is illustrated by the comparison of the two 
reference solutions obtained with the refined refer- 
ence model (solid circles and squares in Figure 9). 
This effect is magnified as the size of the elements 
is increased, which is evident by the results for the 
coarse reference model (solid triangles in the figure). 
The impact of the use of linear strains is decreased 
as the mesh is refined. An enlarged view of the re- 
sults is shown in Figure 9b to delineate better the 
effect of the strain approximation on the transverse 
deflection. 

As was shown in Figure 8, the interface element 
analysis with nonlinear strain approximations (open 
circles in Figure 9) is in excellent agreement with 
the reference solution with nonlinear strain approxi- 
mations. The response obtained with linear strains 
(open triangles in Figure 9) illustrates the stiffening 
of the panel due to the combined effect of the strain 
approximation and the coarse discretization of one of 
the substructures. These results are bounded by the 
coarse and refined reference models using linear 
strain approximations (solid triangles and solid 
squares in the figure). 


6 

American Institute of Aeronautics and Astronautics 





Transverse Deflection, w c 

(b) Enlarged View of Encircled Region 

Figure 9. Effect of Nonlinear Strain 
Approximations on Transverse Deflection. 

Based on these results, nonlinear strain approxi- 
mations are used in the coarse substructure while 
linear strain approximations are used in the refined 
substructure. The results of this analysis with com- 
bined linear and nonlinear strain approximations are 
given by the open squares in Figure 9 and are de- 
noted by coupled (L/NL) in the figure. These results 
are bounded by the solutions using the refined refer- 
ence model with linear strains and nonlinear strains 
(solid squares and circles in the figure, respectively). 
Computations with nonlinear strains are more expen- 
sive than those with linear strains due to the in- 
creased number of operations required in evaluating 
the nonlinear terms in the strain-displacement rela- 
tions. Hence, a trade-off can be made between the 
level of mesh density and the level of nonlinear ap- 
proximation. For this benchmark application, the 
results obtained with linear strains in the most re- 
fined substructure and nonlinear strains in the coarse 
substructure were more accurate than the results ob- 
tained with the refined reference model using linear 


strains only. These results demonstrate the unique 
analysis capability provided by the interface tech- 
nology. Although in this example, the corotational 
formulation provided the nonlinear response charac- 
ter in each substructure, this concept may be ex- 
tended to the application of linear/nonlinear substruc- 
turing in which one or more substructures are identi- 
fied as exhibiting nonlinear behavior while the other 
substructures are assumed to exhibit linear behavior. 
In this case, the tangent stiffness matrices are com- 
puted for the nonlinear substructures and assembled 
along with the linear stiffness matrices for the linear 
substructures. 


Hinged Cylindrical Panel 

An isotropic cylindrical panel which is hinged 
on its two straight edges and free on its two curved 
edges (shown in Figure 10) and which is subjected to 
a concentrated load at its center is used to demon- 
strate the capability of the interface technology on a 
structure exhibiting collapse. This panel configura- 
tion is given in the literature and attributed to Sabir 
and Lock. 21 The panel material properties are 
3.10275 kN/mm 2 (450 ksi) for the Young’s modulus 
and 0.3 for the Poisson’s ratio. The panel radius, R, 
is 2540 mm. (100 in.); the panel half-length, L, is 
254 mm. (10 in.); and the half-opening angle, 9, is 
0.1 radians. Two panel thicknesses, t, were consid- 
ered: a thickness of 12.7 mm. (0.5 in.) and a thick- 
ness of 6.35 mm. (0.25 in.), yielding radius-to- 
thickness ratios, R/t, of 400 and 200, respectively. 
Both the thick and thin panels exhibit a limit point 
and snap-through behavior (See Figure 11) as the 
load is increased, and the panels collapse into an 
inverted configuration. In addition, the thin panel 
exhibits a snap-back behavior (See Figure 12). The 
coupled finite element model used for the analysis of 
both panels is the same as that used for the clamped 
cylindrical panel and is shown in Figure 7. The inter- 
face is located along an axial line at the panel half- 
opening angle, 9. 



Figure 10. Hinged Cylindrical Panel 
Subjected to Concentrated Load. 
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The concentrated load is given as a function of 
the transverse deflection at the center of the thick 
and thin panels (point C in Figure 10) in Figure 11 
and Figure 12, respectively. The response curve for 
the thick panel (see Figure 11) indicates the limit 
point and snap-through of the panel as the load is 
increased. The response curve for the thin panel (see 
Figure 12) indicates the snap-back behavior of the 
panel. The results with the interface element for 
both the thin and thick panels are in excellent 
agreement with the reference solution. 



Figure 11. Central Deflection of Thick Hinged 
Cylindrical Panel. 



Transverse Deflection, w c (mm) 


Figure 12. Central Deflection of Thin Hinged 
Cylindrical Panel. 

Cantilevered Plate in Pure Bending 

An isotropic rectangular plate (shown in Figure 
13) subjected to an end bending moment is used to 
demonstrate the capability of the interface technol- 
ogy for a structure exhibiting very large rotations. 
The plate is clamped at its root, and, to emulate 
beam behavior, symmetry conditions are used on the 
long edges. This configuration is often referred to as 
the elastica problem. The plate material properties 
are 120 psi for the Young’s modulus and 0.0 for the 
Poisson’s ratio. The plate length, L, is 10 in., and 
the width, W, and thickness, t, are both 1 in. The 


coupled finite element model is shown in Figure 14 
with the interface located at the panel midwidth. In 
Figure 15, the normalized moment is given as a 
function of the deflection at the loaded end of the 
plate normalized by the plate length. The results 
indicate the very large rotations which the structure 
undergoes, and the interface element analysis is in 
excellent agreement with the reference solution. The 
interface element performed well even for such large 
rotations. 



S ymmet ry 


Symmetry" 



Clamped 


Figure 13. Cantilevered Plate Subjected to Pure 
Bending. 



Figure 14. Finite Element Models for 
Cantilevered Plate. 



Normalized Moment, M/M max 

Figure 15. Deflection of Cantilevered Plate. 


Composite Panel 

A composite rectangular panel (shown in Figure 
16) subjected to axial compression is used to demon- 
strate the capability of the interface technology in 
the postbuckling regime. This panel is clamped at 
the root and simply- supported by knife-edge supports 
on its long edges to prevent it from buckling as a 
wide column. The panel used is a 24-ply orthotropic 
laminate and is denoted Panel C4 in the experimen- 
tal and analytical results reported by Starnes and 
Rouse. 22 The panel is fabricated from unidirectional 
graphite-fiber tapes preimpregnated with 450K cure 
thermosetting epoxy resin. The ply properties are 
19,000 ksi for the longitudinal Young’s modulus, 
1890 ksi for the transverse Young’s modulus, 930 ksi 
for the in-plane shear modulus, 250 ksi for the trans- 
verse shear modulus, and 0.38 for both Poisson’s ra- 
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tios. The laminate stacking sequence for the panel is 
[±45/ O 2 f± 45 /O 2 f± 45/0/90]^, and the nominal 

ply thickness is 0.00551 in. The panel length, L, is 
20 in., and the width, W, is 6.75 in. The overall 
panel thickness, t, is 0.13 in. The reference model 
has 12 elements along the length and 6 elements 
along the width. The panel was observed in the test 22 
to buckle into two longitudinal half-waves and one 
transverse half-wave which resulted in peak stresses 
along its center. Thus, the region of interest for this 
application is in the vicinity of the center of the 
panel, and a refined model is used in that region. 


The coupled finite element model is shown in 
Figure 17. Two interfaces are used in this analysis 
and are depicted by the gray shaded lines in Figure 
17. The compressive load is given as a function of 
the axial deflection at the loaded edge and maxi- 
mum transverse deflection normalized by the panel 
thickness in Figure 18 and Figure 19, respectively. 
The response curves indicate the softening of the 
panel as the load is increased. The postbuckling 
response exhibits large out-of-plane deflections (over 
three times the thickness, see Figure 19). The inter- 
face element analysis is in excellent agreement with 
the reference solution. 


Clamped 


Simply-Supported 



Simply-Supported 


Figure 16, Composite Panel Subjected to 
Compressive Loading. 


Figure 17. 


\ / 


Interfaces 

Finite Element Models for 
Composite Panel. 



Figure 18. Axial End Deflection of Composite 
Panel. 



Figure 19. Maximum Transverse Deflection for 
Composite Panel. 

Concluding Remarks 

Interface technology for geometrically nonlinear 
analysis of multiple connected subdomains has been 
demonstrated. As in the case of the linear formula- 
tion, the subdomains need not be nodally compati- 
ble, and, thus, the need for potentially complex tran- 
sition modeling is eliminated. A hybrid variational 
formulation was utilized to achieve compatibility, in 
a variational sense, between the independently mod- 
eled substructures. The development of the hybrid 
formulation in the form of an element facilitates the 
implementation of this approach in standard finite 
element software packages. The interface element 
stiffness is computed and assembled along with stan- 
dard finite elements allowing for the use of standard 
finite element matrix assemblers and linear and non- 
linear solution procedures. 

The interface technology described herein has 
been demonstrated on several applications which 
exhibit a variety of response characteristics, namely: 
stiffening, softening, snap-through, snap-back, large 
rotations, and postbuckling. The applications used to 
demonstrate the capabilities of the interface technol- 
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ogy for predicting such responses were a clamped 
square plate loaded with a uniform pressure, a 
clamped cylindrical panel loaded with uniform pres- 
sure, a hinged cylindrical panel with a concentrated 
transverse load, a cantilevered plate subjected to an 
end bending moment, and a composite panel loaded 
in axial compression. The results obtained using the 
interface technology were in excellent agreement 
with the reference solutions for all cases. 

The capability of using different orders of strain 
approximations in the independent substructures has 
been demonstrated with the analysis of the clamped 
cylindrical panel. This unique analysis capability 
allows the analyst to specify the order of the strain- 
displacement approximation to be used in each sub- 
structure. Nonlinear strain approximations are com- 
putationally more expensive than linear strain ap- 
proximations due to the increased number of opera- 
tions required in evaluating the nonlinear terms in 
the strain-displacement relation. Hence, a trade-off 
can be made between the level of mesh density and 
the level of nonlinear approximation. In the case 
studied, the results obtained with linear strains in the 
most refined substructure and nonlinear strains in the 
coarse substructure were more accurate than the re- 
sults obtained with the reference model approxi- 
mated with linear strains only. In addition, as often 
done in nonlinear substructuring, the analyst may use 
the interface technology to isolate a region requiring 
nonlinear strain approximations while using the lin- 
ear strain approximations elsewhere. 

The interface technology described herein pro- 
vides a method for predicting the nonlinear response 
of plate and shell problems as well as concomitant 
detailed stress states. The element has been imple- 
mented within a general-purpose finite element code. 
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